SAIC-00/8030:APPAT-263 
October 1 7, 2000 


'^MODELING ULTRAVIOLET EMISSIONS 
NEAR lO" 

NASW-99013 

Final Report, 

June 1, 1999-May 31, 2000 


NASA CONTRACT: NASW-99013 


Principal investigator: 

Jon a. Linker 

Science Applications International Corporation 
10260 Campus Point Drive 
SAN DIEGO, CA 92121-1578 



SAIC-00/8030:APPAT-263 
October 17, 2000 


MODELING ULTRAVIOLET EMISSIONS NEAR lO 
FINAL REPORT: 

6/1/99 - 5/31/00 
1. Introduction 

In this report, we describe work performed under Contract NASW-99013, 
awarded to Science Applications International Corporation, for the period 6/1/99 to 
5/31/00. During this time period, we have investigated the interaction of lo, 
Jupiter's innermost Galilean satellite, with the lo plasma torus, and the role this 
interaction plays in producing ultraviolet (UV) emissions from neutral oxygen and 
sulfur. 

lo, the innermost of Jupiter's Galilean satellites, plays a unique role in the 
Jovian magnetosphere. Neutral material that escapes from lo is ionized to form the 
lo torus, a dense, heavy-ion plasma that corotates with Jupiter and interacts with lo. 
lo supplies not only the torus, but is a major source of plasma for the entire 
magnetosphere. Ionization and charge-exchange of neutrals near lo strongly 
influences the plasma interaction, and lo’s neutral atmosphere plays an important 
role in the generation of currents that couple lo to Jupiter. There have been no i n 
situ measurements of the neutral density near lo, but remote observations of 
neutrals near lo have been performed for many years [Brown, 1974; Trafton et al. 
1974; Brown et al., 1983 and references therein]. 

Recent observations from the Hubble Space Telescope (HST) have shown 
detailed structure in UV emissions from neutral species near lo. Electron-impact of 
the neutrals by the lo torus plasma is the primary mechanism responsible for 
exciting these emissions. Previously, we have modeled the lo plasma environment 
using three-dimensional magnetohydrodynamic (MHD) simulations (Linker et al. 
1998a), and we have shown that the interaction between lo and the plasma torus 
plays an important role in producing the morphology of the observed emissions 
(Linker et al. 1998bc). In the past year, we have extended these studies to use both 
UV observations and Galileo particle and field measurements to investigate the lo 
interaction. 
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2. Intrinsic Magnetization at lo 

A key question raised by results from the original flyby of the Galileo 
spacecraft past lo (December 7, 1995) is whether or not lo is intrinsically magnetized. 
The question is not likely to be settled conclusively with the data from this flyby. 
The first of three new flybys of lo by Galileo was planned for October 11, 1999 
(referred to as 124). In anticipation of this flyby, as well as a coordinated HST 
observing campaign, we used our MHD model to perform further studies of the 
plasma interaction at lo. We considered both magnetized and conducting models 
of lo to see if differences in the predicted UV emission, together with possible 
differences in the predicted plasma and magnetic field signatures for 124, might 
allow us to distinguish more clearly between the models. 

In the original flyby of lo (referred to as JO), the Galileo spacecraft flew through 
lo's wake with a closest approach distance of 898 km. Both conducting and 
magnetized models of lo could qualitatively match many of the features observed in 
the Galileo data, including the large depression in magnetic field magnitude that 
was observed [Linker et ah, 1998]. The 124 trajectory passed by lo upstream and on 
the flanks of lo. From the point of view of simple potential (current-free) magnetic 
field models, one might expect that for an Ionian dipole moment approximately 
aligned with the local (Jovian) magnetic field one would see a depression in the 
local magnetic field both well upstream and downstream of lo. This would be in 
contrast to what is expected for a conducting lo, where draping of the magnetic field 
upstream of lo causes the local magnetic field strength to increase while the currents 
driven in the conducting region reduce the field strength downstream of lo. 
However, the situation at lo is more complicated, because even if lo does have an 
intrinsic magnetic field, significant ion pickup must also be occurring and this will 
cause slowing of the flow and draping of the magnetic field. It turns out that for 
reasonable ion pickup rates, the region of depressed magnetic field upstream of lo in 
the magnetized case is closely confined near lo and not likely to be penetrated 
significantly by the 124 flyby. Figure 1 shows a summary of predictions for the 
plasma and magnetic field measurements during the 124 flyby that were presented at 
the October 1999 DPS meeting [Linker et al., 1999]. Given the poor constraints on 
parameters of the model (e.g., ion pickup rate, extent of the collisional ionosphere, 
strength and tilt of any intrinsic field) our results suggest that conclusive proof or 
refutation of an lo intrinsic magnetic field are not likely from the Galileo 
measurements during this encounter. 
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Comparison of Conducting and Magnetized Models of lo 
Predictions for the Galileo 124 Flyby 
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FigureL Comparsion of predictions for the 124 flyby from conducting and magnetized models of lo. 
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To investigate whether UV observations from HST might reveal significant 
differences between the magnetized and conducting models of lo, we used the 
plasma properties from the MHD computations to simulate the observed emission 
for different viewing geometries. For collisional excitation of neutrals by electron 
impact, the apparent emission rate Ccin be written as an integral over the line of 
site: 

* (') “f' 

where 


e (/) 




exp 



and Tin is the density of the neutral emitter, K is a constant related to the collision 
strength, is the electron density, Tg is the electron temperature (in eV), and AE is 
the energy of the emitted photon (in eV). In the MHD computations, the plasma 
properties near lo are computed self-consistently for a given neutral density 
distribution (n„ ). To compute the emission, we use nn (assumed to be spherically 
symmetric), the plasma density, and the plasma temperature from the MHD 
simulation. Te is assumed to be a fixed fraction of the total plasma temperature. 
We orient the simulation domain such that the plasma flow direction is along the 
direction of corotation and the background magnetic field is aligned with the local 
direction of Jupiter's magnetic field (at lo) at the time of the observation. 

Because the UV emissions are a line-of-sight integrated quantity, their 
properties can depend strongly on the viewing geometry. Figure 2 shows the 
definition of lo longitudes with respect to views from Earth, and Figure 3 shows 
how the how the simulated emission appears for different lo phases for both a 
conducting model of lo (top panels) and a magnetized model of lo (bottom panels). 
When lo is near longitude 0 (behind Jupiter) or 180 (in front of Jupiter), both models 
predict that the brightest emission occurs on the downstream side of lo. Near 
elongation, the models predict that the brightest emission occurs on the flanks of lo; 
this pattern of emission was observed in images taken with the Space Telescope 
Imaging Spectrograph (STIS). 

While there are systematic differences between the emission brightness 
predicted by the two models, the emission morphology is very similar. Because the 
brightness can be very sensitive to model assumptions (background plasma 
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Figure 3. The variation of simulated planeofsky UV emission at 1356 A (Ol) as a function of lo’s 
longitudefor a conducting (top 4 panels) and magnetized (bottom 4 panels.. The lo torus plasma corotates 
with Jupiter at 57 km/s faster than lo rotates. Therefore, at longitude 0, lo’s wake is to the left and upstream 
of lo is to the right. At longitude 90, lo’s wake is visible. At longitude 180, lo’s wake is to the right and 
upstream of lo is to the left. At longitude 270, lo’s upstream side is visible. The magnetized and conducting 
models show similar morphologies, but different emission brightnesses. 
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temperature and neutral density, for example), it will be difficult to distinguish 
between conducting and magnetized models on the basis of UV emission alone. 

3. Comparison of Simulated and Calibrated UV Images 

The techniques we have developed for computing simulated images can be 
used to compare our results with specific observations. Figure 4 shows an example 
of such a comparison for an HST STIS observation of OI (1356A) taken on October 
14, 1997. Figure 4a shows the local magnetic field geometry for the time of the 
observation. We use the viewing geometry at the time of the observation and 
compute the emission as in Figure 3. Figure 4b and 4c show a comparison between 
the STIS image and a simulated image generated with the plasma parameters from 
the MHD model. In the STIS image, we see that there is more intense emission on 
lo's flanks, there is a tilt of the emission in the plane (lo north is up), and there is 
higher emission on the northern part of lo relative to the south. These properties 
are all reproduced in the simulated image. 

HST STIS Observation: October 14, 1997 
lo CML = 269, lo System III = 84, Jupiter CML = 173 
Jupiter 




Figure 4, Comparison of a simulated emission image with a specific HST observation, (a) The geometry 
of a STIS observation taken on October 14, 1997. The orientation of a Jovian magnetic field line passing 
near lo is shown, (b) An HST STIS observation at 1356 A (OI). The brightest emission is shown in red, 
the faintest emission is blue, (c) A simulated planeofsky image of UV emission at 1 356 A (OI) for the 
geometry shown. The image was created using the plasma density and temperature predicted by the 
MHD simulation, and a spherically symmetric neutral density. 
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We have performed a number of comparisons of our simulated images 
developed using the MHD model with HST observations (Linker et al. 1998bc, 
Linker, 1999; McGrath and Linker 1999). Many of the morphological features of the 
observations are reproduced in the simulated images, indicating the importance of 
lo's plasma environment in producing the UV emissions from neutral species. 
However, these comparisons were qualitative in nature, in that the observations 
were not all calibrated and only a relative brightness was computed for the 
simulated images. To perform a more quantitative comparison, we have developed 
a calibrated image for a specific STIS observation and also computed the absolute 
brightness. 

Figure 5 shows a comparison of a calibrated HST STIS image with the 
absolute brightness predicted by the MHD simulation. Note that the model matches 
the morphological features of the observation reasonably well. In particular, the 
high intensity "spots" on the limbs and their tilt relative to lo's equator are present 
in both the model and the observations, and the spot on the Jupiter side of lo has 
slightly less intensity. When the images are scaled exactly the same, we see that the 
peak emission predicted by the model is about 50% less than that obtained in the 
calibrated image. Given the uncertainties in many parameters (such as the neutral 
density at lo), the agreement is quite good. 

The results from our MHD computations of the lo interaction and 
comparisons with both Galileo data and HST observations of UV emission show 
that the models can reproduce many features of the observations. Our models 
suggest that much of the structure in the UV emissions from neutrals species near 
lo can be explained as a production of the interaction of the lo torus plasma with lo's 
atmosphere and exosphere. Even with an assumption of a spherically symmetric 
exosphere, the "spot-like" features on the flanks of lo that have been consistently 
observed when lo is near elongation are reproduced. While lo's atmosphere may in 
actuality have significant departures from spherical symmetry, our results suggest 
that particular asymmetric features of the atmosphere are not essential for 
producing the observed emission features. 
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HST STIS Observation: October 14, 1997 

lo CML = 243, lo System III = 357, Jupiter CML = 60 
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Figure 5. Comparison of a calibrated HST STIS image with absolute brightness predicted from the MHD 
model (a) The geometry of the STIS observation. The orientation of a Jovian magnetic field line 
passing near lo is shown. Note that lo was not yet at elongation (in contrast to the case shown in Figure 
4). (b) An HST STIS observation at 1356 A (OI). The image on the left shows the approximate position 
of lo as a red circle, (c) A simulated planeofsky image of UV emission at 1356 A (OI) for the 
geometry shown. The absolute brigthness was computed using the plasma density and temperature 
predicted by the MHD simulation. The image on the left is scale betwen 140 and 1600 Rayleighs (the 
approximate range for the predicted emissiion, while the image on the right is scaled exactly the same as 
the STIS image (300 to 2000 Rayleighs). 
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